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according  to  strain  increments  conjugate  to  Cauchy  stress  is  found  to  be  the  most  stringent 
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order  constants  tend  to  result  in  more  unstable  behavior,  regardless  of  the  sign  of  the  pres¬ 
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much  more  intrinsically  stable  than  second-  and  third-order  elastic  models  when  Poisson’s 
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1.  Introduction 

Stability  of  elastic  solids  under  finite  deformation  has 
been  the  subject  of  numerous  studies,  with  early  work  on 
intrinsic  stability  of  crystals  due  to  Born  (1940).  Nonlinear 
elastic  anisotropic  solids  (e.g.,  single  crystals)  have  been 
analyzed  in  a  number  of  works  (Hill,  1 975 ;  Hill  and  Milstein, 
1977;  Milstein  and  Hill,  1979;  Wang  et  ah,  1993,  1995; 
Morris  and  Krenn,  2000),  as  have  nonlinear  elastic  isotropic 
solids  (Hill,  1957;  Rivlin,  1974;  Rivlin  and  Beatty,  2003). 

Various  criteria  for  stability  of  elastic  solids  have  been 
proposed  in  the  literature,  beginning  with  work  of  Born 
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(1940)  who  associated  stability  with  a  positive  definite 
stiffness  measure  and  local  convexity  of  internal  energy 
expressed  in  terms  of  a  Lagrangian  Green  strain  measure. 
In  a  real  physical  system,  the  appropriate  choice  of  stability 
criterion  depends  on  the  method  of  static  incremental  load 
application  (e.g.,  dead  loading  in  one  or  more  directions 
(Rivlin,  1974))  and  any  constraints  associated  with  bound¬ 
ary  conditions,  and  such  a  criterion  may  not  correspond  to 
Born’s.  When  considering  “intrinsic”  or  “internal’’  stability 
of  a  unit  cell  or  unit  cube  of  a  given  material,  from  the  con¬ 
tinuum  viewpoint  or  using  atomic  theory,  the  proper 
choice  of  stability  measure  is  ambiguous  when  the  precise 
loading  mechanism  is  left  unspecified.  Different  choices  of 
conjugate  stress-strain  measures  (i.e.,  different  general¬ 
ized  coordinates  and  conjugate  forces)  can  lead  to  different 
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local  convexity  conditions  and  different  intrinsic  stability 
criteria  (Hill,  1975;  Hill  and  Milstein,  1977).  As  stated  by 
Milstein  and  Hill  (1979),  consistency  of  classical,  unique, 
and  environment-dependent  stability  criteria  with  intrin¬ 
sic  stability  or  convexity  arguments  such  as  Born’s  requires 
a  special  environment  in  which  loads  can  be  varied  to  fol¬ 
low  the  material  during  any  disturbances  while  fixing  val¬ 
ues  of  conjugate  forces,  which  in  turn  are  non-unique  since 
they  depend  on  the  choice  of  generalized  coordinates  (i.e., 
the  choice  of  strain  measure).  This  load  environment  used 
to  probe  stability  (e.g.,  via  virtual  deformations  at  a  fixed 
stress)  need  not  correspond  to  the  loading  program  used 
to  achieve  the  stressed  equilibrium  state  from  which  sta¬ 
bility  is  tested  (Hill  and  Milstein,  1977;  Milstein  and  Hill, 
1979). 

Intrinsic  stability  requirements  considered  in  the  pres¬ 
ent  work  lead  to  restrictions  on  various  incremental  tan¬ 
gent  elastic  moduli;  stability  measures  of  this  sort 
applied  to  superposed  deformation  variations  in  any  direc¬ 
tion  may  require  that  an  energy  function  be  locally  convex, 
or  at  least  that  a  particular  stiffness  matrix  (i.e.,  Hessian)  be 
positive  definite,  at  the  current  equilibrium  state.  For  non¬ 
linear  materials  of  high  symmetry  (e.g.,  isotropic  solids) 
under  simple  deformation  paths,  such  intrinsic  stability 
requirements  can  often  be  stated  succinctly  in  terms  of 
restrictions  on  a  few  strain-dependent  elastic  coefficients. 
In  linear  elastic  solids,  such  requirements  degenerate  to 
the  usual  constitutive  constraints  of  positive  bulk  and 
shear  moduli.  As  noted  above  and  explained  in  Hill  and 
Milstein  (1977)  and  Parry  (1978),  intrinsic  stability  mea¬ 
sures  like  those  considered  herein  do  not  depend  on  the 
load  environment,  but  do  depend  on  the  choice  of  conju¬ 
gate  stress-strain  measures.  Comparisons  among  intrinsic 
criteria  involving  different  Lagrangian  strain  measures 
were  derived  for  generic  isotropic  elastic  solids  under 
positive  principal  stresses  (Parry,  1978). 

The  distinction  between  material  instability  and  struc¬ 
tural  instability  should  be  noted.  As  defined  in  the  present 
work,  material  instability  correlates  with  intrinsic  instabil¬ 
ity,  and  depends  only  on  material  properties  and  loading 
protocol.  Intrinsic  stability  criteria  are  local  since  they  con¬ 
sider  only  homogeneous  stress/deformation  states  up  to 
the  onset  of  instability.  In  nonlinear  materials  (e.g.,  nonlin¬ 
ear  elastic,  elastic-plastic,  or  damaged  solids),  the  onset  of 
material  instability  depends  on  strain,  but  in  linear  elastic 
solids,  material  stability  is  independent  of  strain  and 
simply  requires  positive  definiteness  of  the  tensor  of  elas¬ 
tic  constants.  In  contrast,  structural  instability  criteria  are 
global  rather  than  local,  depending  on  geometry  of  the 
body.  The  stress/strain  state  may  be  inhomogeneous  prior 
to  onset  of  structural  instability.  Structural  instability  may 
occur  even  if  the  material  is  intrinsically  stable,  and  can  be 
induced  by  loads  even  in  linear  elastic  materials.  A 
representative  example  is  buckling  of  a  slender  column 
under  compression. 

Intrinsic  stability  properties  of  solids  under  large  com¬ 
pressive  stress  or  pressure  are  of  interest  for  applications 
in  ballistics,  impact  phenomena,  and  earth  and  planetary 
sciences.  Elastic  instability  may  signal  the  onset  of  failure 
or  localization  phenomena,  e.g.,  slip,  fracture,  or  phase 
transformations  (Hill,  1975;  Wang  et  al.,  1993;  Morris 


and  Krenn,  2000).  Hard  materials  of  lower  symmetry  such 
as  quartz  (Gregoryanz  et  al.,  2000),  silicon  carbide 
(Clayton,  2010),  and  boron  carbide  (Clayton,  2012)  exhibit 
a  decrease  in  certain  shear  elastic  stiffness  components 
with  increasing  pressure.  At  high  pressures,  this  tendency 
may  lead  to  the  onset  of  instability  and  subsequent  amor- 
phization  (Chen  et  al.,  2003),  which  in  the  case  of  ceramic 
materials  hinders  performance  in  ballistic  applications.  In 
contrast  to  these  polyatomic  ceramic  materials,  most 
elemental  engineering  materials  demonstrate  increasing 
shear  moduli  with  increasing  pressure  (Guinan  and 
Steinberg,  1974),  which  would  tend  to  enhance  rather  than 
diminish  internal  stability  at  large  compressions. 

Nonlinear  elastic  models  of  anisotropic  single  crystals 
(Wallace,  1972;  Teodosiu,  1982;  Clayton,  2011)  typically 
assume  a  strain  energy  function  written  as  a  Taylor  poly¬ 
nomial  in  Green  (Lagrangian)  elastic  strain 
E  =  i(FTF-l),  with  F  the  deformation  gradient.  Such  a 
model,  when  terms  of  up  to  third  order  are  maintained 
(i.e.,  second-  and  third-order  elastic  constants)  provides 
reasonably  accurate  descriptions  of  stresses  and  wave 
propagation  for  moderate  compressions  (Thurston,  1974; 
Clayton,  2009);  however,  Eulerian  strain  measures  may 
be  more  accurate  for  extreme  pressures  (Weaver,  1976; 
Jeanloz,  1989). 

The  present  work  focuses  primarily  on  isotropic  elastic 
solids  of  third  order.  Such  materials  are  described  by  two 
independent  second-order  elastic  constants  and  three 
independent  third-order  elastic  constants  (Murnaghan, 
1937;  Teodosiu,  1982).  Because  of  the  limited  number  of 
constants,  systematic  study  of  effects  of  choices  of  con¬ 
stants  on  intrinsic  stability  for  simple  monotonic  deforma¬ 
tion  paths  is  tractable,  and  is  undertaken  in  this  work. 
Third-order  constants  are  related  explicitly  to  pressure 
derivatives  of  bulk  and  shear  moduli  in  the  reference  state 
(Thurston  et  al.,  1966;  Guinan  and  Steinberg,  1974; 
Teodosiu,  1982);  experimental  data  (Guinan  and  Steinberg, 
1974;  Steinberg,  1982)  thus  provide  realistic  bounds  on 
combinations  of  third-order  constants.  Choices  of  third-or¬ 
der  constants  yielding  a  decreasing  shear  stiffness  with 
increasing  compression  provide  insight  into  behavior  of 
aforementioned  ceramic  materials  demonstrating  shear 
instabilities  at  high  pressure  (Gregoryanz  et  al.,  2000;  Chen 
et  al.,  2003).  It  is  noted  that  single  crystals  of  such  materi¬ 
als  are  highly  anisotropic  (e.g.,  trigonal  symmetry:  six  sec¬ 
ond-order  and  fourteen  third-order  elastic  constants)  and 
systematic  study  of  effects  of  varying  all  elastic  constants 
individually  on  stability  is  intractable.  For  comparison, 
intrinsic  stability  of  a  class  of  compressible  neo-Hookean 
solids  (Simo  and  Pister,  1984),  which  demonstrate  a 
strongly  increasing  bulk  modulus  with  pressure,  is  also 
considered. 

This  paper  is  organized  as  follows.  Requisite  quantities 
associated  with  internal  stability  are  derived  in  Section  2. 
Intrinsic  stability  of  third-order  elastic  solids  in  terms  of 
three  different  criteria  from  the  literature  (Born,  1940;  Hill, 
1975;  Wang  et  al.,  1993)  is  analyzed  in  Section  3.  For  each 
criterion,  minimum  eigenvalues  of  a  particular  tangent 
stiffness  matrix  are  examined  for  different  choices  of  sec¬ 
ond-  and/or  third-order  elastic  constants  for  an  element 
of  material  undergoing  spherical  deformation,  uniaxial 
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strain,  and  simple  shear.  Implicitly  assumed  in  this  analy¬ 
sis  are  availability  of  so-called  "passive”  loading  environ¬ 
ments  (Hill  and  Milstein,  1977;  Milstein  and  Hill,  1979) 
that  make  possible  application  of  any  of  the  three  criteria 
mentioned  already,  permitting  unstable  modes  in  any 
strain  direction  that  may  or  may  not  be  coaxial  with  the 
imposed  primary  deformation  path  up  to  the  checkpoint 
for  instability.  Complementary  analysis  of  such  internal 
stability  of  neo-Hookean  solids  is  also  undertaken. 
Although  some  general  constraints  have  been  derived  else¬ 
where  for  cubic  crystals  under  pressure  (Wang  et  al.,  1993) 
and  neo-Hookean  solids  under  dead  loading  (Rivlin,  1974; 
Rivlin  and  Beatty,  2003),  effects  of  ranges  of  elastic  con¬ 
stants  on  all  of  the  considered  intrinsic  stability  criteria 
for  either  of  these  two  kinds  of  material  models  have  not, 
to  the  authors’  knowledge,  been  reported  elsewhere. 

Notational  conventions  of  continuum  mechanics  are 
used.  Vectors  and  higher-order  tensors  are  written  in  bold 
italics,  scalars  and  scalar  components  in  italics.  Stiffness 
tensors  and  matrices  are  written  in  serif  fonts  (e.g.,  A.  B, 
etc.)  All  are  referred  to  a  fixed  Cartesian  frame,  with  refer¬ 
ential  components  in  capitals  and  spatial  components  in 
lower-case.  Greek  indices  are  used  for  Voigt  notation. 
Summation  applies  over  repeated  indices. 


Stresses  are 


2  F 


dW 

9C 


EdW 

Fw 


=  FS. 


(2.8) 


Second-order  tangent  elastic  moduli  are 


A(F) 


dP 

OF 


d2W 
OFdF ’ 


C  (E) 


as 

dE 


a2w  _  o2w 

OEdE  ~  ^  dCDC ' 


(2.9) 


These  are  related  by  Simo  and  Pister  (1984)  and  Clayton 
(2011)  as 

A  iju  =  FuFkKCi]ia  +  SjL5ik  (2.10) 

and  possess  the  symmetry  properties 

f^ijkL  =  Afeljj,  CljKL  =  CklIJ  =  CjlKL  =  Qjlk-  (2.1 1) 


Three  intrinsic  or  internal  stability  criteria,  all  fre¬ 
quently  encountered  in  various  solid  mechanics,  physics, 
or  materials  science  literature,  are  addressed  in  what  fol¬ 
lows.  First,  consider  stability  with  respect  to  allowable 
variations,  from  a  possibly  stressed  equilibrium  state,  in 
Green  strain,  SE.  The  following  internal  stability  criterion 
(also  denoting  local  convexity)  is  usually  attributed  to  Born 
(1940): 

SS  ■  SE  =  SE  ■  C  ■  SE  =  SExCvpSEp  >  0,  (2.12) 


2.  Nonlinear  elasticity  and  intrinsic  stability  measures 

2.3.  General  theory 

Let  x  and  X  denote  spatial  and  reference  coordinates. 
The  deformation  gradient  is  the  invertible  two-point 
tensor 

F  =  dx/dX,  Fij  =  dXi/dXj,  (2.1 ) 


where  [C^]  is  the  symmetric  6x6  matrix  corresponding  to 
C  in  Voigt  notation  (a ,p=  1,2 _ 6).  The  actual  net  poten¬ 

tial  energy  change  (internal  work  to  second  order  minus 
virtual  work  at  fixed  S)  is  half  the  quantity  on  the  left  side 
of  inequality  (2.12)  in  the  context  of  a  trapezoidal  rule  for 
quadrature  (Hill,  1975);  inclusion  or  omission  of  a  factor  of 
l  does  not  change  the  stability  criterion.  A  necessary  and 
sufficient  condition  for  (2.12)  to  hold  under  arbitrary  non¬ 
zero  SE  is 


where  i.J  =  1,2,3.  By  the  polar  decomposition  theorem, 


F  =  RU  R'=Rt  U  =  Ut,  (2.2) 

where  U  is  positive  definite.  The  Jacobian  determinant 
giving  the  ratio  of  current  to  initial  volume  is 

J  =  detF  =  detU  >  0.  (2.3) 

The  Green  (Lagrangian)  strain  tensor  and  deformation 
tensor  are 

E  =  i(C-l),  Es  =  l(Ca-da)-,  (2.4) 

C  =  FtF  =  U2 ,  Clj=Fk,Fkj  =  UIKUjK.  (2.5) 


Cauchy  stress  a  =  <tt,  first  Piola-Kirchhoff  stress  P,  and 
second  Piola-Kirchhoff  stress  (i.e.,  thermodynamic  ten¬ 
sion)  S  =  SJ  are  related  by 

a  -J  'PF1  -J  1 FSF1  av  =J  'PlKFiK 
=J  'F„SLKF]K.  (2.6) 


det[Ca/j]  >  0.  (2.13) 

At  instability,  [Ca/(]  has  null  or  negative  eigenvalue(s),  with 
corresponding  eigenmodes  <3Ea.  When  the  current  configu¬ 
ration  is  taken  as  reference, 

Cijki(F)  =J  ^FuFjjFitKFuCijKL,  detj  =  F Kf  F ^  SErl  (2.14) 

and  since  J  >  0,  criterion  (2.12)  and  (2.13)  is  equivalent  to 

SEijCijklSEkl  =  JSeijCijkiSeid  >  0  x=>  det[ca?]  >  0.  (2.15) 

Nonzero  eigenvalues  of  c  and  C  may  be  different,  but  the 
critical  deformation  F*  at  which  either  first  exhibits  a  zero 
eigenvalue  is  the  same.  Criteria  (2.12)  and  (2.15)  are  equiv¬ 
alent  to  Eqs.  (4.1 )— (4.3)  of  Hill  (1975).  Although  early  cal¬ 
culations  (Born,  1940;  Misra,  1940)  focused  on 
unstressed  lattices,  this  criteria  has  often  been  posited  for 
stressed  materials. 

A  second  internal  stability  criterion  can  be  expressed  as 
the  inequality  (Hill,  1957,  1975) 

<5P  SF  =  SF  •  A  •  <5F  >  0.  (2.16) 


Cauchy  pressure  is  p  =  —  i  cri{. 

For  a  hyperelastic  material  with  properties  independent 
of  X,  strain  energy  per  unit  reference  volume  W  is  of  the 
general  form 

W  =  W[C(F)]  =  W[E(F)]  =  W(F).  (2.7) 


where  SF  is  an  admissible  variation  in  deformation  gradi¬ 
ent.  Criterion  (2.16)  is  equivalent  to  Eq.  4.20  of  Hill 
(1975)  applied  in  the  context  of  dead  loading.  For  cases 
in  which  variations  are  restricted  to  strain  only, 
i5F  — >  SU  =  (SU)T,  and  (2.16)  reduces  to,  considering  all 
possible  strain  increments, 
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SU  ■  A  •  SU  =  SUaA^SU!!  >  0  det[Aa/)]  >  0,  (2.17) 

where  [A^]  is  the  symmetric  6x6  matrix  corresponding  to 

Aym  =  ^  (AijiddiiditK  +  AjikidjjdkK  +  AjjiKdiidu  +  Ajmdjjdn). 

(2.18) 

When  the  current  configuration  is  taken  as  reference, 

a,j id(F)  =  7  1  FjyF/LA,yfct,  <5e,j  =  ^(FK-SFjK  +  FK)SFiK)  (2.19) 

and  since  J  >  0,  criterion  (2.17)  can  also  be  expressed  as 

SUijAjklSUkl  =Jdeljaijkideki  >  0  det[a^]  >  0.  (2.20) 

As  noted  in  Clayton  (2012),  this  is  equivalent  to  the  so- 
called  GCN+  condition,  listed  as  constitutive  inequality 
(8A.21)  of  Wang  and  Truesdell  (1973). 

A  third  intrinsic  stability  criterion  is  (Wang  et  al„  1993; 
Morris  and  Krenn,  2000) 

So  ■  Se  =  Se  ■  B  •  Se  >  0,  (2.21) 

where  Se  is  a  nonzero  strain  increment  defined  by 

<5  e  =  F-TSEF-\  de,j  =  FK^dEKLF,/  (2.22) 

and  work  conjugate  to  Cauchy  stress  since 

3W  =  S  ■  SE  =  P  ■  SF  =Jo  ■  Se.  (2.23) 

Incremental  tangent  modulus  giving  change  in  Cauchy 
stress  with  respect  to  strain  is  (Wallace,  1972;  Wang 
et  al„  1993;  Morris  and  Krenn,  2000;  Clayton,  2012) 

B,jw  =  Cjjid  +  ^(0'iilSji  +  Gudjk  +  ajkdu  +  ffjiSjh  -  CTijdki  -  Ukldij). 

(2.24) 

Tensor  B  in  (2.24)  exhibits  full  Voigt  symmetry: 

Byw  =  B  kiij  —  Bji^/  =  Bjik.  (2.25) 

Correspondingly,  denoting  [B^]  the  symmetric  6x6  form 
of  B,  internal  stability  holds  with  respect  to  Cauchy  stress 
if  and  only  if 

det[B^]  >  0.  (2.26) 

Condition  (2.21)  is  identical  to  Eq.  (30)  of  Morris  and  Krenn 
(2000),  and  the  corresponding  instability  criterion  is  iden¬ 
tical  to  Eq.  (2.26)  of  Wang  et  al.  (1995).  Symmetrized 
tangent  stiffness  B  is  identical  in  the  current  work  and 
these  other  works  (Wang  et  al.,  1993;  Wang  et  al.,  1995; 
Morris  and  Krenn,  2000;  Clayton,  2012);  notation  used  in 
Wang  et  al.  (1993)  and  Wang  et  al.  (1995)  for  “Lagrangian” 
strain  corresponds  to  an  increment  in  Lagrangian  strain 
measured  with  respect  to  a  stressed  reference  state  and 
is  identical  to  (2.22)  of  the  present  work,  as  derived  explic¬ 
itly  by  Wallace  (1967).  A  criterion  similar  to  (2.21)  was 
declared  to  be  in  accord  with  the  “exact”  condition  for  clas¬ 
sical  stability  of  hydrostatically  loaded  cubic  crystals  with 
moduli  dictated  by  pairwise  atomic  interactions  (Milstein 
and  Hill,  1979). 

Although  the  first  equality  in  (2.21)  does  not  strictly 
apply  for  cases  in  which  rotation  of  the  material  1, 
condition  (2.26)  has  been  interpreted  (Morris  and  Krenn, 
2000)  as  a  necessary,  but  not  always  sufficient,  intrinsic 


stability  criterion  for  arbitrary  (i.e„  unconstrained)  load 
increments.  Utility  of  criteria  incorporating  variations  in 
rotation  such  as  (2.16)  with  SF  ¥=  (SF)T  as  conditions  for 
the  onset  of  structural  changes  in  materials  is  question¬ 
able,  since  even  the  zero  stress  state  is  a  boundary  point 
of  the  (un)stable  domain  when  rotations  arise  (Hill, 
1975).  Noting  that  under  rotations  of  the  spatial  frame, 
Se  transforms  objectively  (similarly  to  the  symmetric 
deformation  rate  tensor)  but  Str  generally  does  not,  the 
leftmost  side  of  (2.21)  is  not  necessarily  objective  with  re¬ 
spect  to  changes  among  rotating  spatial  frames  of  refer¬ 
ence.  This  deficiency  in  the  B  stability  criterion  has 
apparently  not  been  emphasized  in  prior  works  that  em¬ 
ploy  it  (Wang  et  al.,  1993;  Wang  et  al„  1995;  Morris  and 
Krenn,  2000),  and  no  attempt  is  made  to  address  it  further 
here.  When  Str  =  -Sp  1  is  a  pressure  increment,  objectivity 
of  (2.21 )  does  hold.  Furthermore,  term  Se  ■  B  ■  Se  following 
the  equality  in  (2.21),  and  hence  criterion  (2.26),  is  invari¬ 
ant  under  change  of  spatial  basis  since  B  transforms  like  a 
fourth-order  tensor. 

Minimum  eigenvalues  of  symmetric  6x6  matrix  forms 
of  C,  A,  and  B  are  labeled  Ac,  Aa,  and  AB.  Since  in  the  unde¬ 
formed  and  stress  free  reference  state  C  =  A  =  B,  the  min¬ 
imum  eigenvalue  of  all  of  these  matrices  in  said  reference 
state  is  denoted  A0. 

Considered  later  in  Section  3  are  homogeneous  defor¬ 
mations  of  a  volume  element  (e.g.,  a  unit  cube)  of  uniform 
material.  Each  of  the  C,  A,  and  B  intrinsic  stability  criteria  is 
treated  as  equally  plausible  since  the  particular  means  by 
which  incremental  traction  is  applied  to  such  a  volume  un¬ 
der  virtual  deformations  from  the  homogeneously  de¬ 
formed  state  is  left  unspecified.  For  a  specific  situation 
involving  dead  loading,  for  which  traction  per  unit  refer¬ 
ence  area  and  hence  P  is  held  fixed  during  stretch  variation 
U,  the  A  criterion  (2.17)  would  be  most  appropriate.  When 
the  loading  mechanism  fixes  the  Cauchy  stress,  e.g.,  as  in 
conventional  hydrostatic  compression  or  tension  (Milstein 
and  Hill,  1979;  Wang  et  al.,  1993;  Wang  et  al.,  1995),  then 
B  criterion  (2.21 )  would  be  most  appropriate,  as  discussed 
in  Milstein  and  Hill  (1979)  and  Morris  and  Krenn  (2000). 
The  C  criterion  (2.12)  does  not  correlate  readily  with  any 
real  load  control  mechanism  (since  there  is  no  traction  vec¬ 
tor  associated  with  S  that  is  independent  of  deformation), 
but  is  considered  in  this  paper  because  of  its  historical  sig¬ 
nificance  (Born,  1940;  Misra,  1940;  Hill,  1975)  and  wide¬ 
spread  use  in  physics  and  materials  science  (Gregoryanz 
et  al.,  2000).  Viewed  another  way,  a  comparison  of  convex¬ 
ity  criteria  as  generalized  coordinates  (i.e.,  strain  mea¬ 
sures)  vary  is  analogous  to  a  comparison  of  classical 
elastic  stability  criteria  as  the  load  environment  varies 
(Parry,  1978). 

Consideration  of  finite  strain  and  correct  work  conju¬ 
gate  stress-strain  measures  is  essential  for  meaningful 
study  of  deformation-induced  material  instability.  These 
different  measures  enter  the  C,  A,  and  B  criteria  defined 
above.  In  contrast,  in  linear  elasticity  theory,  different 
stress  tensors  (S.P.a)  and  strain  tensors  ( E,U,e )  become 
indistinguishable,  and  intrinsic  stability  does  not  depend 
on  strain  since  tangent  stiffness  is  constant;  e.g.,  for  an  iso¬ 
tropic  linear  elastic  solid  the  local  stability  requirement 
degenerates  to  (3.5). 
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2.2.  Third-order  elasticity 


Strain  energy  density  written  as  a  Taylor  series  in  Green 
strain  E  is 

W  =  W  +  C  yEy  +  —  CijKLEijEkl  +  CijklmneijeklEmn  +  •  •  • 


=  W°+  — 
dE,j 


F  1  &u 

o  U  2dEgdE„. 


EI]EKL 


1 

+  = 


6  dE,]dEKLdEMN 


E=0 

EijEklEmn  +  • 


(2.27) 


Vanishing  energy  and  stress  in  the  undeformed  reference 
state  leads  to  W°  =  0  and  Cj)  =  0,  and  henceforth  restrict¬ 
ing  attention  to  third-order  elasticity, 

W  =  2  E-ijklEijEkl  +  g  (~ijklmneUEklemn 

=  2CvA£/'  +  6C^£/1Ez.  (2.28) 

Stress  S  and  tangent  modulus  C  are 

Sa  =  C^E,  + 1  C^E*.  C*  =  C°„  +  C  °pxEx.  (2.29) 

Second-  and  third-order  constants  in  Voigt  notation  are  C°p 
and  C°flr  For  homogeneous  isotropic  materials,  elastic 
constant  tensors  are  of  the  form  Teodosiu  (1982) 

Cj/m.  =  +  G^ikSjl  +  &  ilSjk),  (2.30) 


G°jklmn  =  ^[SijSklSmn]  +  P[Sij{SKmSln  +  SknSlm) 

+  Skl(SimSjn  +  5in8]m)  +  Smn(&i k^jl  +  <5jl<S/k)] 

+  X\ Sik(SjmSln  +  SjNSLm)  +  <>jl(8imSkn  +  <5/ nSkm) 

+  5/l(<5jm<5/cn  +  SjnSkm)  +  Sjk(SimSln  +  <5;n<5im)]- 

(2.31) 

The  theory  of  third-order  elastic  isotropic  solids  with  five 
independent  elastic  constants  was  developed  by  Murna- 
ghan  (1937),  Second-order  elastic  constants  obey  the 
familiar  relations 

7  =  C?2,  G=C^4  =  i(C?1-C?2),  K 

-J  +  SC'  ”-2i <lW  <232> 

where  C.K,  and  v  are  shear  modulus,  bulk  modulus,  and 
Poisson’s  ratio  in  the  reference  state.  Third-order  elastic 
constants  obey 

a=C?23’  P  =  C?44  =  2  (Cll2  _  C°23),  Z  =  C456 


=  8(C?„-3C?12  +  2C' 


0  1 
123/  * 


(2.33) 


Second-order  compliance  components  =  [C^]  1  are 
^~0  |  ^-0 

M?i=- 


C°  4-  C° 
Ml  +  M2 


(C?,-C°2)(C°1+2C°2)! 


M"  = 


— ru 

M: 


(C?,-C?2)(C?1+2C?2)’ 


it  follows  that  m“4  =  2(Mjj  -  Mj2)  =  l/C^  =  1/G.  Deriva¬ 
tives  of  second-order  moduli  Ca/(  with  respect  to  pressure 
are,  at  the  undeformed  reference  state  (i.e.,  F  =  1  =^>  £  =  0), 


,  _  dCxp 
dP  f-i 


(2.35) 


where  in  full  tensor  notation  (Thurston  et  al.,  1966), 


c>  — 

^  I  Ilf  I  — 


dQjKL  0EMn  dSpQ 


®Kl  9Emn  dSpQ  dp 


=  -C 


o 

IJKLMNIV,MNPP 


Mr 


(2.36) 


For  third-order  isotropic  materials, 

Cn  =  —  3^(^ni  2C112) , 

C'2  =  -  JR  (2C?12  +  C?23).  (2.37) 

leading  to  pressure  derivatives  of  thermodynamic  bulk  and 
shear  moduli  at  the  reference  state: 


K'  =  -g^(C?„+6C?12+2C?23),  C' 

=  -^(C?„-C?23).  (2.38) 

Conditions  for  which  stiffness  tensor  C  remains  isotropic 
under  finite  deformation  have  been  recently  discussed 
(Fuller  and  Brannon,  2011).  Full  Cauchy  symmetry 
(Clayton,  2011)  corresponds  to  X  =  C  and  oc  =  p  =  %. 


2.3.  Neo-Hookean  elasticity 


For  a  compressible  neo-Flookean  material,  strain  energy 
density  is 

W  =  W°  +  ^C(trC-3)+/(/) 

=  W°  + 1  G(F,/Fi7  -  3)  +/(/),  (2.39) 

where  function  /accounts  for  compressibility.  For  the  par¬ 
ticular  kind  of  neo-Hookean  solid  considered  here,  which 
demonstrates  vanishing  stress  and  energy  in  the  reference 
state  (Sirno  and  Pister,  1984), 

W  =  lc(trC-3)  -GlnJ  +  ^l(\nJ)2.  (2.40) 

Noting  that  8\n]/dC  =  1C-1,  stress  S  and  tangent  modulus 
C  are 

Sl,  =  2^=Cd,]  +  (;.\nJ-C)Cl]\  (2.41) 

acij 


d2W 

CnflCvi 


(2.42) 


In  the  reference  state,  C1  =1,  7  =  1,  and  (2.30)  and  (2.42) 
are  consistent;  relations  among  second-order  constants  for 
neo-Hookean  elasticity  obey  (  2.32).  Considering  spherical 
deformation  of  the  form  C/;  =  J2/3dy  such  that  cr,j  =  -pSy, 
pressure  derivatives  of  bulk  and  shear  moduli  at  the  refer¬ 
ence  state  are 


(2.34) 
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!('  = 


a  = 


K 


4 

+  3 


and  analyzed  for  respective  Green  elastic  and  compressible 
neo-Hookean  solids  in  Section  3.2  and  Section  3.3. 

3.2.  Green  elasticity 


Extension  of  this  neo-Hookean  model  to  elastic-plastic 
materials  and  expressed  in  general  curvilinear  coordinates 
can  be  found  elsewhere  (Clayton,  2012). 

3.  Intrinsic  stability  under  particular  deformation  paths 

Homogeneous  deformations  of  an  element  of  isotropic 
material  are  considered.  Minimum  eigenvalue  ratios 
Ac/Ao,Aa/Ao,  and  AB/A0  are  compared  for  fully  imposed 
spherical,  uniaxial,  and  shearing  paths,  recalling  that  a  zero 
eigenvalue  corresponds  to  the  onset  of  intrinsic  instability, 
e.g.,  loss  of  local  convexity  associated  with  a  given  stiffness 
tensor  for  the  current  deformed  state.  Eigenvalues  and  cor¬ 
responding  orthonormal  eigenvectors  (i.e.,  eigenmodes) 
are  computed  numerically  as  the  deformation  is  incre¬ 
mented  along  each  path.  Eigenmodes  are  not  restricted  to 
align  with  the  primary  (imposed)  deformation  path;  when 
critical  modes  do  not,  then  possibility  of  bifurcation  from 
the  primary  path  is  implied. 

At  the  conclusion  of  each  applied  increment  of  defor¬ 
mation,  it  is  assumed  that  a  (possibly  fictitious)  loading 
apparatus  is  available  enabling  virtual  variations  in  strain 
(SE,  SU,  or  Se)  in  any  direction  while  holding  the  conjugate 
stress  measure  (S.P,  or  a,  respectively)  fixed  (Hill,  1975; 
Hill  and  Milstein,  1977).  Such  conditions  enable  formal 
application  of  the  intrinsic  stability  criteria  defined  in 
Section  2.1.  A  similar  approach  was  taken  in  Clayton 
(2012)  in  a  study  of  internal  stability  of  a  particular 
material  of  lower  symmetry  (boron  carbide  ceramic), 
where  following  (Wang  et  al.,  1995),  internal  stability 
was  referred  to  as  “mechanical  stability”. 

3.1.  Imposed  deformation  paths 

Prescribed  deformation  is  parameterized  as  follows  in 
terms  of  a  single  scalar  £,  letting  {ei,e2,e3}  denote  an 
orthonormal  basis.  Spherical  (and  hence  hydrostatic  stress 
states  for  isotropic  materials): 

F  =  U  =  (  1+£)1/31,  J=l  +  £;  (3.1) 

Uniaxial  strain  along  e ^ : 

F  =  U  =  l  +  eei®ei,  J=  l+£;  (3.2) 

Simple  shear  in  the  et  -  e2  plane: 

F  =  RU  =  1  +  ee,  ®e2,  J=  1;  (3.3) 


r  2 
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r  2 
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(4+£2)1/2 

(4+e2)1/2 
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(4+e2)1/2 

(4+£2),/2 
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l  o 

0 

1_ 

0 

0 

lj 

(3.4) 

Effects  of  various  values  of  elastic  constants  on  internal 
stability  under  the  above  deformation  paths  are  computed 


Considered  first  are  materials  for  which  all  third-order 
elastic  constants  vanish,  i.e.,  St.  Venant-Kirchhoff  solids. 
In  this  case,  Qjkl  =  C°JKL  =  constant,  and  internal  stability 
constraints  in  the  unstressed  and  undeformed  reference 
state  are 

det[Ca/!]  >  0  4=^  A0  =  Ac  >  0  4=>  I(  >  0  and  G  >  0. 

(3.5) 

These  can  also  be  associated  with  local  strain  energy  con¬ 
vexity  at  the  reference  state.  The  rightmost  conditions  in 
(3.5)  imply  -1  <  v  <  \.  Noting  that  G  =  i(A  -  1),  at  a  given 
value  of  v,  energy  W,  stresses,  and  tangent  elastic  coeffi¬ 
cients  are  all  proportional  to  X.  Hence  Aa/A0  and  AB/A0  de¬ 
pend  only  on  v  and  e  for  a  particular  deformation  path. 
Ratio  Ac/Ao  —  1  regardless  of  v  and  e  and  is  of  no  interest. 
Considered  in  Fig.  1  are  results  for  -1  <  v  <  J.  Intrinsic 
instability  occurs  as  A/A0  — >  0,  i.e.,  as  a  curve  intersects 
the  horizontal  axis  of  each  figure. 

As  shown  in  Fig.  1(a)  and  (b)  for  spherical  deformation, 
as  v  increases,  stability  in  tension  (e  >  0)  increases,  and  in 
compression  (e  <  0)  decreases.  The  B  stability  criterion  of 
(2.26)  tends  to  be  more  stringent  than  the  A  criterion  of 
(2.17),  i.e.,  a  zero  eigenvalue  is  more  often  attained  at  a 
smaller  magnitude  of  £  for  the  former,  particularly  for 
v  >  0.  As  shown  in  Fig.  1(c)  and  1(d)  for  uniaxial  strain, 
if  v  ^  0  then  instability  may  occur  in  both  tension  and 
compression.  As  shown  in  Fig.  1(e)  and  (f)  for  simple  shear, 
stability  increases  with  increasing  v,  with  Aa/A0  and 
Ab/A0  remaining  positive  for  v  ^  0  and  £  ^  0.8.  The  B  cri¬ 
terion  provides  a  more  stringent  measure  of  stability  than 
the  A  criterion  for  shear  deformation.  As  discussed  more  in 
the  Appendix,  as  the  material  tends  towards  incompress¬ 
ibility  (v  — >  0.5),  intrinsic  stability  declines  rapidly  under 
spherical  or  uniaxial  compression. 

Even  though  third-order  constants  are  zero  for  results 
shown  in  Fig.  1,  the  theory  is  still  geometrically  nonlinear 
and  mathematically  valid  at  large  magnitudes  of  finite 
strain  parameter  £,  i.e.,  results  shown  do  not  correspond 
to  classical  linear  elasticity  theory,  which  does  not  distin¬ 
guish  among  Piola-Kirchhoff  and  Cauchy  stresses,  for 
example.  Vanishing  third-order  constants  corresponds  to 
K'  =  G'  =  0  in  (2.38),  which  as  will  be  noted  later,  is  an 
exceptional  case  that  does  not  hold  for  most  elemental 
polycrystalline  solids.  Results  in  Fig.  1  are  important,  how¬ 
ever,  because  St.  Venant-Kirchhoff  theory  is  widely  used, 
especially  in  the  context  of  finite  elastic-plastic  theories 
of  metal  plasticity  (Rashid  and  Nemat-Nasser,  1992; 
Clayton  and  McDowell,  2004;  Clayton,  2006).  Results  also 
provide  a  basis  of  comparison  for  calculations  reported 
next  wherein  stiffness  C  is  not  constant. 

Considered  next  are  cases  for  which  one  or  more  third- 
order  constants  (i.e.,  a,  / 1 ,  %  or  C°n,  C°12,  C°23)  may  be  non¬ 
zero.  Attention  is  focused  on  ranges  of  values  physically 
realistic  for  polycrystalline  engineering  materials  (e.g., 
metals,  ceramics,  and  minerals).  Ultrasonic  measurements 
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Fig.  1.  Normalized  minimum  eigenvalues  of  stiffness  tensors  A  and  B  for  isotropic  second-order  elasticity  (CjL,  =  0):  (a)  spherical  compression,  A  (b) 
spherical  compression,  B  (c)  uniaxial  compression,  A  (d)  uniaxial  compression,  B  (e)  shear,  A  (f)  shear,  B. 


and  various  high  pressure  experiments  suggest  that  the 
pressure  derivative  of  the  bulk  modulus  in  the  reference 
state  obeys  2  <  K1  <7  for  most  such  materials,  with 
3  <  K'  <  6  typical  (Steinberg,  1982;  Jeanloz,  1989).  Experi¬ 
mental  data  also  suggest  that  the  pressure  derivative  of  the 
thermodynamic  shear  modulus  in  the  reference  state 
1  <  G'  <  4  for  pure  elements  ((Guinan  and  Steinberg, 
1974),  correcting  for  G'  differing  from  the  definition  in  that 
reference  by  G/3I<  +  1 ).  Also  of  high  interest  in  the  context 
of  instability  at  high  pressure  are  polyatomic  materials 


with  G'  <  0  such  as  silicon  carbide  (C44  ss  -0.2)  (Clayton, 
2010)  and  boron  carbide  (-5  <  C44  <  -2)  (Clayton, 
2012).  Means  of  constructing  effective  isotropic  third- 
order  constants  for  polycrystals  from  anisotropic  single 
crystal  constants  are  summarized  in  Mason  and  Maudlin 
(1999)  and  references  therein. 

In  this  work,  third-order  constants  probing  the  domain 
0  </<:'<  10  and  -IO^G'^15  are  considered,  as  listed  in 
Table  1.  In  the  third-order  elastic  analysis,  attention  is 
restricted  to  v  =  \  (i.e„  2  =  G  =  |/C),  a  sensible  value 
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Table  1 

Values  of  third-order  elastic  constants  investigated  for  intrinsic  stability. 
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representative  of  many  metals,  ceramics,  and  minerals.  As 
discussed  by  Gilman  (2003),  v  =  i  (i.e.,  C°2  =  C^)  is  charac¬ 
teristic  of  the  alkali  metals  and  ionic  crystals  such  as  alkali 
halides  that  tend  towards  Cauchy  symmetry.  Covalently 
bonded  materials,  e.g.,  many  ceramics,  tend  to  have 
while  engineering  metals  tend  to  have  slightly  larger  v, 
e.g.,  the  mode  is  closer  to  i  for  face  centered  cubic  metals. 
Pure  solid  elements  span  the  range  0.02  <  v  <  0.45 
(Gilman,  2003).  Some  values  for  particular  materials  are 
now  quoted.  Sodium  chloride  (alkali  halide,  rocksalt  struc¬ 
ture)  and  alumina  (ceramic,  trigonal  structure)  have 
v  fa  0.24  (Clayton,  2009).  Two  covalently  bonded  materials 
are  silicon  carbide  (v  fa  0.16,  hexagonal  structure  (Clayton, 
2010))  and  diamond  (v  fa  0.07,  cubic  structure  (McSkimin 
et  al.,  1972)).  Most  single  crystals  are  anisotropic;  values 
of  v  mentioned  thus  far  correspond  to  isotropic  polycrys¬ 
tals.  Two  metals  whose  single  crystals  have  low  anisotropy 
are  tungsten  (v  fa  0.28,  body  centered  cubic  structure 
(Clayton,  2005))  and  magnesium  (v  f»  0.27,  hexagonal 
structure  (Clayton  and  Knap,  2011)).  Concrete,  essentially 
a  porous  mixture  of  minerals,  has  via  0.15  (Clayton, 
2008),  though  variations  with  composition  are  probable. 
In  summary,  v  =  |  is  deemed  a  reasonable  compromise  of 
values  mentioned  above. 

At  fixed  v,  it  can  be  shown  that  energy  W,  stresses,  and 
tangent  elastic  coefficients  depend  only  on  X  and  C°^/X, 
and  A0  depends  only  on  X  (linearly).  Hence  Ac/A0,  Aa/A0, 
and  Ab/A0  depend  only  on  C°ft/2  and  e  for  a  particular 
deformation  path.  For  example,  for  Case  1  with  v  =  l, 


Case  1.  [Fig.  2(a),  K'  >  0  and  G' >  0]  The  material 
becomes  increasingly  stable  (unstable)  in  com¬ 
pression  (tension)  with  increasing  K'  and  G'. 
The  C,  A,  and  B  criteria  give  similar  trends,  with 
B  more  stringent  in  compression  and  C  more 
stringent  in  tension. 

Case  2.  [Fig.  2(b),  K'  >  0  and  G'  =  0]  The  material 
becomes  increasingly  unstable  in  compression 
with  increasing  /<".  A  change  in  critical  eigen¬ 
mode  is  observed  in  tension  for  K’  =  8.3,  leading 
to  instability  at  e  >  0.12.  The  B  criterion  tends  to 
be  more  stringent,  especially  in  compression. 

Case  3.  [Fig.  2(c),  K'  >  0  and  G'  <  0]  The  material 
becomes  increasingly  unstable  in  compression 
with  decreasing  G'.  The  B  criterion  tends  to  be 
more  stringent,  especially  in  compression. 

Case  4.  [Fig.  2(d),  K'  >  0  and  G'  =  0]  Results  are  identical 
to  those  for  Case  2.  This  demonstrates  the  fact 
that  for  spherical  loading  of  isotropic  third- 
order  elastic  materials,  internal  stability 
depends  only  on  K,G,IC,  and  G',  and  not  on  all 
three  independent  third-order  elastic  constants 
considered  individually. 

Case  5.  [Fig.  2(e),  IC  =  0  and  G'  >  0  or  G'  <  0)  The  mate¬ 
rial  becomes  increasingly  stable  in  compression 
(tension)  with  positive  (negative)  G'.  The  C,  A, 
and  B  criteria  give  similar  trends,  with  B  more 
stringent  in  compression  and  C  more  stringent 
in  tension. 
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where  0  >  C°u/X  >  -150  in  Table  1. 

Shown  in  Fig.  2  are  results  for  spherical  deformation;  in 
this  deformation  program  stress  is  hydrostatic,  i.e., 
a  =  -pl.  Each  subfigure  shows  minimum  normalized 
eigenvalues  of  C,  A,  and  B  for  two  different  sets  of  third- 
order  elastic  constants  within  prescribed  domains  for  each 
of  Cases  1-5  listed  in  Table  1.  Recall  that  A/A0  =  1  at  e  =  0 
for  any  stiffness  matrix,  and  that  A/A0  =  0  at  the  onset  of 
intrinsic  instability  as  defined  in  the  present  work. 
Noteworthy  observations  are  as  follows: 


Results  in  Fig.  2  confirm  a  particular  example  of  a  more 
general  comparison  theorem  derived  in  Parry  (1979)  [Eq. 
(36)  of  that  reference]:  for  a  fixed  material  model,  convex¬ 
ity  with  respect  to  stretch  (i.e.,  A  stability)  is  weaker  than 
convexity  with  respect  to  Green  strain  (i.e.,  C  stability) 
when  Sn  =  S22  =  S33  =  -_/1/3p  3=  0. 

Shown  in  Fig.  3  are  results  for  uniaxial  compression/ 
tension;  in  this  deformation  program  stress  usually  in¬ 
cludes  significant  hydrostatic  and  deviatoric  components. 
Noteworthy  observations  are  as  follows: 

Case  1.  [Fig.  3(a),  IC  >  0  and  G' >  0]  The  material 

becomes  more  unstable  in  tension  with  increas¬ 
ing  K'  and  G'.  The  C,  A,  and  B  criteria  give  similar 
results  for  the  critical  strain  at  instability. 

Case  2.  [Fig.  3(b),  K'  >  0  and  G'  =  0]  The  material 

becomes  increasingly  unstable  in  compression 
and  tension  with  increasing  K'.  The  B  criterion 
tends  to  be  more  stringent,  especially  in 
compression. 
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(a)  (b) 


Fig.  2.  Normalized  minimum  eigenvalues  of  stiffness  tensors  C.  A,  and  B  for  isotropic  third-order  elasticity  and  spherical  deformation:  (a)  case  1,  Cjn  /  0 
(b)  case  2,  C°12  t6  0  (c)  case  3,  C°23  ^  0  (d)  case  4,  =  C°12  =  C?23  ¥•  0  (e)  case  5,  C?,,  =  -2C°23  ^  0. 


Case  3.  [Fig.  3(c),  K'  >  0  and  G'  <  0]  The  material 
becomes  increasingly  unstable  in  compression 
and  tension  with  decreasing  G'.  The  domain  of 
intrinsic  stability  is  small  for  G'  =  -6.3,  i.e. 
-0.03  <  £  <  0.05.  The  C,  A,  and  B  criteria  give 
similar  results  in  compression. 

Case  4.  [Fig.  3(d),  K'  >  0  and  G'  =  0]  Unlike  what  was 
observed  for  spherical  deformation,  results  for 
uniaxial  strain  differ  from  those  in  Case  2. 


However,  like  Case  2,  the  material  becomes 
increasingly  unstable  in  compression  and  ten¬ 
sion  with  increasing  /<". 

Case  5.  [Fig.  3(e),  K' =  0  and  G' >  0  or  G'  <  0]  The 
domain  of  intrinsic  stability  shifts  towards  com¬ 
pression  (tension)  with  positive  (negative)  G'. 
The  C,  A,  and  B  criteria  give  similar  trends,  with 
B  more  stringent  in  compression  and  C  more 
stringent  in  tension. 
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Fig.  3.  Normalized  minimum  eigenvalues  of  stiffness  tensors  C.  A.  and  B  for  isotropic  third-order  elasticity  and  uniaxial  strain:  (a)  case  1,  C /  0  (b)  case 
2,  C?12  ^  0  (c)  case  3,  C?23  s4  0  (d)  case  4,  =  C?12  =  C?23  z4  0  (e)  case  5,  =  -2C?23  ^  0. 


Shown  in  Fig.  4  are  results  for  simple  shear;  in  this 
deformation  program  stress  includes  comparatively 
large  deviatoric  components,  though  hydrostatic  stress 
does  not  generally  vanish  due  to  geometric  nonlinear¬ 
ity,  i.e„  Kelvin  and  Poynting  effects  (Eringen,  1962), 
and  may  be  large  as  a  result  of  large  strain  or  large 
material  nonlinearity.  Noteworthy  observations  are  as 
follows: 

Case  1.  [Fig.  4(a),  K' >  0  and  G'  >  0]  The  material 
becomes  more  unstable  at  large  shear  with 


increasing  K'  and  G' .  The  C,  A,  and  B  criteria  give 
similar  results  for  the  critical  strain  at 
instability. 

Case  2.  [Fig.  4(b),  K'  >  0  and  G'  =  0]  The  material 
becomes  increasingly  unstable  in  shear  with 
increasing  /<".  The  B  criterion  is  the  most 
stringent. 

Case  3.  [Fig.  4(c),  I<'  >  0  and  G'  <  0]  The  material 
becomes  increasingly  unstable  in  shear  with 
decreasing  G'.  The  C,  A,  and  B  criteria  give  simi¬ 
lar  results,  with  C  most  stringent. 
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Fig.  4.  Normalized  minimum  eigenvalues  of  stiffness  tensors  C,  A,  and  B  for  isotropic  third-order  elasticity  and  shear  deformation:  (a)  case  1,  C°u  /-  0  (b) 
case  2,  C?12  0  (c)  case  3,  Cj23  sz  0  (d)  case  4,  =  C°12  =  C°23  ?z  0  (e)  case  5,  C211  =  -2C°23  ?z  o. 


Case  4.  [Fig.  4d,  K'  >  0  and  C'  =  0]  Unlike  what  was 
observed  for  spherical  deformation,  results  for 
simple  shear  differ  from  those  in  Case  2.  The 
material  becomes  increasingly  unstable  with 
increasing  /<",  but  very  large  strains  ( e  >  0.45) 
are  attained  prior  to  instability. 

Case  5.  [Fig.  4(e),  K'  =  0  and  C'  >  0  or  G'  <  0]  Compari¬ 
son  with  Fig.  1(e)  and  1(f)  demonstrates  that 
the  material  becomes  increasingly  unstable 
with  increasing  magnitude  of  G',  regardless  of 
sign  of  G'.  The  case  with  G'  <  0  has  a  slightly 


larger  critical  strain  (e*  ss  0.12)  than  the  case  with 
G'  >  0  (£*  ==  0.11).  The  C,  A,  and  B  criteria  give 
similar  predictions  for  the  minimum  eigenvalue. 

By  inspection,  Ac  does  not  depend  on  rotation  R  enter¬ 
ing  (3.4).  For  results  presented  in  this  paper,  it  was  also 
found  that  critical  strain  e*  at  instability  corresponding  to 
Aa  =  0  or  AB  =  0  did  not  change  if  stretch  U  was  pre¬ 
scribed  rather  than  RU  in  (3.4);  i.e.,  rigid  rotation  did  not 
affect  loss  of  internal  stability  in  simple  shear.  Specifically, 
Ab  does  not  depend  at  all  on  rotation  R  since  B  is  an 
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objective  fourth  order  tensor,  and  while  Aa  can  depend  on 
R,  such  dependence  was  found  to  be  inconsequential  at  ap¬ 
plied  shear  strains  up  to  the  onset  of  internal  instability. 

3.3.  Neo-Hookean  elasticity 

In  the  reference  state,  CIJKl  =  C°JKL  and  internal  stability 
conditions  are  again  the  familiar  constraints 

det[C°,]  >  0  A0  >  0  I<  >  0  and  G  >  0.  (3.6) 


The  rightmost  conditions  in  (3.6)  again  imply  -1  <  v  <  J. 
Noting  that  G  =  1(A  -  1),  at  a  given  value  of  v,  it  follows 
that  strain  energy  W,  stresses,  and  tangent  elastic  coeffi¬ 
cients  are  all  proportional  to  X.  Hence  Ac/A0,  Aa/A0,  and 
Ab/A0  depend  only  on  v  and  e  for  a  particular  prescribed 
deformation  path.  Considered  in  Fig.  5  are  results  for 
-1  <  v  <  \.  These  results  are  analogous  to  those  of  Fig.  1, 
but  with  a  compressible  neo-Hookean  model  rather  than 
a  second-order  Green  elasticity  model. 


(b) 


Fig.  5.  Normalized  minimum  eigenvalues  of  stiffness  tensors  A  and  B  for  neo-Hookean  elasticity:  (a)  spherical  compression,  A  (b)  spherical  compression,  B 
(c)  uniaxial  compression,  A  (d)  uniaxial  compression,  B  (e)  shear,  A  (f)  shear,  B. 
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As  shown  in  Fig.  5(a)  and  (b)  for  spherical  deformation, 
as  v  increases,  stability  in  tension  (e  >  0)  decreases,  and  in 
compression  (e  <  0)  increases.  The  A  and  B  stability  criteria 
of  (2.17)  and  (2.26)  yield  notably  different  results.  For  non¬ 
negative  v,  comparison  with  Fig.  1(a)  and  (b)  shows  that 
the  neo-Hookean  material  tends  to  be  much  more  stable 
in  compression  than  the  second-order  elastic  solid.  As 
shown  in  Fig.  5(c)  and  (d)  for  uniaxial  strain,  for 
-0.25  $  v  ^  0.25,  the  material  remains  intrinsically  stable 
over  the  large  domain  of  strain  studied  (-0.5  ^  e  ^  0.5), 
though  the  A  [Fig.  5(c)]  and  B  [Fig.  5(d)]  criteria  provide 
different  results.  As  shown  in  Fig.  5(e)  and  (f)  for  simple 
shear,  Aa  and  AB  remain  positive  for  v  5=  -0.5  and 
£  ^  1.0.  Trends  for  Ac  are  qualitatively  similar  to  those 
for  Aa  and  AB  and  are  shown  in  Fig.  6  for  v  = 

3.4.  Summary  and  discussion  of  critical  eigenmodes 

Critical  strain  £*  and  corresponding  eigenmode  compo¬ 
nents  5e*  at  the  onset  of  B  instability  (AB  =  0  at  £*)  are 
listed  in  Table  2  for  representative  cases  considered  in 
Section  3.3.  Results  are  given  for  v  =  \.  The  B  stability  crite¬ 
rion  is  presumably  of  highest  interest  since  results  in 
Section  3.3  demonstrated  it  to  be  equally  or  more  stringent 
than  C  or  A  criteria  for  compressive  and  shear  loading 
paths,  for  most  values  of  elastic  constants  considered. 

For  spherical  or  uniaxial  loading,  (+)  denotes  tension 
with  £*  >  0,  and  (-)  denotes  compression  with  e*  <  0. 
Absent  rows  denote  that  a  zero  eigenvalue  is  not  attained 
for  that  case.  Extra  rows  for  a  particular  case  denote  a 
duplicate  null  eigenvalue  and  its  corresponding  critical 
eigenmode,  one  for  each  multiplicity  of  zero  eigenvalue. 
Recall  that  an  eigenmode  with  =  5e2  =  ^  #  0  and 
5e\  =  Sel  =  <5eB  =  0  corresponds  to  a  volumetric  instability; 
conversely,  an  eigenmode  with  &e\  +  Se2  +  de2  =  0  corre¬ 
sponds  to  a  shear  instability.  An  intrinsic  volumetric  insta¬ 
bility  is  associated  with  a  degenerate  incremental  bulk 
modulus,  and  in  a  real  material  might  signify  cavitation 
(Wang  et  al.,  1993;  Wang  et  al„  1995)  or  mode  I  fracture. 
An  intrinsic  shear  instability  is  associated  with  a  degener¬ 
ate  incremental  shear  modulus,  and  in  a  real  material 
might  signify  localized  slip,  twinning  (Rosakis  and  Tsai, 


1994),  or  mode  11/111  fracture.  An  eigenmode  that  is  neither 
purely  volumetric  nor  pure  shear  is  referred  to  as  a  mixed 
instability.  The  following  observations  from  Table  2 
regarding  internal  instabilities  are  noteworthy: 

Second-order  elasticity.  Instability  under  spherical  com¬ 
pression  corresponds  to  five  simultaneous  shearing  modes. 
Instabilities  under  uniaxial  compression  and  simple  shear 
are  mixed. 

Third-order  elasticity,  Case  1.  Instability  under  spherical 
extension  corresponds  to  five  simultaneous  shearing 
modes.  Instabilities  under  uniaxial  extension  and  simple 
shear  are  mixed. 

Third-order  elasticity,  Case  2.  Instability  under  spherical 
extension  is  volumetric.  Instability  under  spherical  com¬ 
pression  corresponds  to  five  simultaneous  shearing  modes. 
Instability  under  uniaxial  extension  is  a  shearing  mode. 
Instabilities  under  simple  shear  and  uniaxial  compression 
are  mixed. 

Third-order  elasticity,  Case  3.  Instability  under  spherical 
extension  is  volumetric.  Instability  under  spherical  com¬ 
pression  corresponds  to  five  simultaneous  shearing  modes. 
Instability  under  uniaxial  extension  is  mixed.  Instability 
uniaxial  compression  corresponds  to  two  simultaneous 
shearing  modes.  Instability  under  simple  shear  is  mixed. 

Third-order  elasticity,  Case  4.  Instability  under  spherical 
extension  is  volumetric.  Instability  under  spherical  com¬ 
pression  corresponds  to  five  simultaneous  shearing  modes. 
Instability  under  uniaxial  extension,  uniaxial  compression, 
and  simple  shear  are  mixed. 

Third-order  elasticity,  Case  5.  Instability  under  spherical 
compression  corresponds  to  five  simultaneous  shearing 
modes.  Instabilities  under  uniaxial  extension,  uniaxial 
compression,  and  simple  shear  are  mixed. 

Neo-Hookean  elasticity.  No  intrinsic  instabilities  are  ob¬ 
served  or  even  approached  for  spherical  loading 
( |e|  ^  0.8),  uniaxial  loading  (|e|  ^  0.8),  or  simple  shear 
loading  (e  ^  1 )  for  v  =  i;  therefore,  |£*|  »  0.8. 

Cases  2,  3,  and  4  demostrate  volumetric  intrinsic  stabil¬ 
ities  in  hydrostatic  extension  which  are  coaxial  with  the 
prescribed  primary  deformation  path.  All  other  instabili¬ 
ties  occur  in  directions  that  differ  from  the  direction  of  pri¬ 
mary  loading,  meaning  virtual  strain  increments  are  not 


e 

Fig.  6.  Normalized  minimum  eigenvalues  of  stiffness  tensors  C  and  B  for  neo-Hookean  elasticity  with  v  ). 
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Table  2 

Critical  e*  (positive  in  tension  or  shear)  and  corresponding  eigenmodes  Se)  for  B  instability,  v  =  j. 


Model  (Case) 

K' 

G' 

Loading 

e* 

Se- 

Se  \ 

Se\ 

< 

Se-5 

Se- 

Second-order 

0 

0 

spherical  (-) 

-0.40 

0.21 

-0.79 

0.58 

0 

0 

0 

-0.79 

0.21 

0.58 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

uniaxial  (-) 

-0.31 

0.96 

-0.19 

-0.19 

0 

0 

0 

shear 

1.2 

0.27 

0.68 

0.047 

0 

0 

-0.68 

Third-order  (1)  (C?n/A  = 

-150) 

10 

15 

spherical  (+) 

0.042 

0.21 

-0.79 

0.58 

0 

0 

0 

-0.79 

0.21 

0.58 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

uniaxial  (+) 

0.017 

0.94 

-0.24 

-0.24 

0 

0 

0 

shear 

0.034 

0.33 

0.32 

-0.22 

0 

0 

0.86 

Third-order  (2)  (C°12/A  = 

-25) 

10 

0 

spherical  (+) 

0.10 

0.58 

0.58 

0.58 

0 

0 

0 

spherical  (-) 

-0.22 

0.41 

0.41 

-0.82 

0 

0 

0 

0.71 

-0.71 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

uniaxial  (+) 

0.077 

0 

0.71 

-0.71 

0 

0 

0 

uniaxial  (-) 

-0.067 

-0.85 

0.37 

0.37 

0 

0 

0 

shear 

0.088 

-0.25 

-0.28 

0.46 

0 

0 

0.80 

Third-order  (3)  (C?23//l  = 

-75) 

10 

15 

~T 

spherical  (+) 

0.10 

0.58 

0.58 

0.58 

0 

0 

0 

spherical  (-) 

-0.067 

-0.27 

-0.53 

0.80 

0 

0 

0 

-0.77 

0.62 

0.15 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

uniaxial  (+) 

0.045 

-0.35 

0.66 

0.66 

0 

0 

0 

uniaxial  (-) 

-0.026 

0 

0.71 

-0.71 

0 

0 

0 

0 

0 

0 

1 

0 

0 

shear 

0.044 

0.11 

0.092 

-0.53 

0 

0 

0.84 

Third-order  (4)  (C?n/A  = 

-50/3) 

10 

0 

spherical  (+) 

0.10 

0.58 

0.58 

0.58 

0 

0 

0 

spherical  (-) 

-0.22 

0.71 

-0.71 

0 

0 

0 

0 

0.41 

0.41 

-0.82 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

uniaxial  (+) 

0.094 

0.50 

0.61 

0.61 

0 

0 

0 

uniaxial  (-) 

-0.21 

0.92 

-0.28 

-0.28 

0 

0 

0 

shear 

0.43 

0.48 

0.60 

0.37 

0 

0 

-0.52 

Third-order  (5)  (C?n/A  = 

+200/3) 

0 

-10 

spherical  (-) 

-0.054 

-0.81 

0.47 

0.34 

0 

0 

0 

0.078 

0.67 

-0.74 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

uniaxial  (+) 

0.11 

-0.053 

0.71 

0.71 

0 

0 

0 

uniaxial  (-) 

-0.038 

0.97 

-0.18 

-0.18 

0 

0 

0 

shear 

0.080 

0.26 

0.29 

0.22 

0 

0 

-0.89 

Neo-Hookean 

26 

15 

7 

5 

any 

»0.8 

aligned  with  prescribed  deformations  of  Section  3.1.  In  a 
real  material,  a  local  strain  perturbation  due  to  a  defect 
in  the  material  might  be  sufficient  to  trigger  a  localization 
or  failure  mode  in  such  a  direction,  presuming  any  con¬ 
straints  imposed  by  the  environment  permit  such 
phenomena. 

Apparent  trends  in  intrinsic  stability  of  isotropic  solids 
under  primary  deformation  paths  of  spherical  compres¬ 
sion,  uniaxial  strain  compression,  and  simple  shear  have 
been  deduced  in  this  work.  These  are  perhaps  the  simplest 
fundamental  deformation-controlled  loading  paths,  and 


correlate  with  real  experiments  such  as  diamond  anvil  cell 
compression,  shock  compression,  and  direct  shear  loading, 
respectively.  Other  fundamental  loading  protocols  not 
studied  here  include  uniaxial  stress  (Morris  and  Krenn, 
2000)  and  all-around  dead  loading  (Rivlin,  1974),  which 
have  been  analyzed  elsewhere,  but  not  for  third-order 
Green  elastic  and  compressible  neo-Hookean  models  stud¬ 
ied  in  the  present  work.  Results  cannot,  in  general,  be 
immediately  translated  to  other  loading  paths  possible  in 
generic  applications.  As  demonstrated  by  Parry  (1980), 
no  universally  weakest  stability  criterion  exists.  For  a  fixed 
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Table  3 

Eigenvalue  ratios,  B  stability,  spherical  and  uniaxial  compression,  v  >  0.4. 


Spherical,  e  =  -0.05 

Spherical,  e  =  -0.10 

Uniaxial,  e  =  -0.05 

Uniaxial,  e  =  -0.10 

(Ab/A0)v_0.4 

0.7437 

0.4737 

0.7084 

0.4324 

rl 

-4.2079 

-18.0217 

-4.2040 

-18.7845 

r2 

15.7529 

12.1098 

16.4888 

12.6963 

Tf 

10.3652 

10.1742 

10.4046 

10.2124 

rB 

10.0352 

10.0171 

10.0389 

10.0208 

r5 

10.0035 

10.0017 

10.0039 

10.0021 

rB 

10.0004 

10.0002 

10.0004 

10.0002 

r? 

10.0000 

10.0000 

10.0000 

10.0000 

material,  the  most  stringent  criterion  for  stability  depends 
on  the  loading  mode.  For  a  fixed  loading  mode,  the  most 
stringent  criterion  depends  on  material  properties.  Com¬ 
parison  theorems  among  different  criteria  that  depend  on 
different  Lagrangian  coordinates  have  been  derived  (Parry, 
1978;  Parry,  1979);  as  noted  already  in  the  text  describing 
Fig.  2,  one  such  theorem  agrees  with  the  present  analysis  of 
C  and  A  criteria  for  third-order  elastic  bodies  under  hydro¬ 
static  tension. 

4.  Conclusions 

Intrinsic  stability  of  several  kinds  of  nonlinear  elastic 
isotropic  solids  has  been  analyzed  for  homogeneous  defor¬ 
mation  paths,  specifically  spherical  extension-compres¬ 
sion,  uniaxial  tension-compression,  and  simple  shear. 
Loss  of  internal  stability  according  to  the  C,  A,  or  B  criterion 
corresponds  to  first  appearance  of  a  null  eigenvalue  of  an 
incremental  stiffness  matrix  associated  with  second  Pio- 
la-Kirchhoff,  first  Piola-Kirchhoff,  or  Cauchy  stress,  respec¬ 
tively.  Effects  of  choices  of  second-  and  third-order  elastic 
constants  on  stability  have  systematically  studied  for  sec¬ 
ond-  and  third-order  Green  elastic  solids  and  compressible 
neo-Hookean  solids.  For  most  cases  investigated  here, 
especially  those  involving  compression,  intrinsic  stability 
according  to  increments  in  Cauchy  stress  (i.e.,  B  criterion) 
is  found  to  be  the  most  stringent.  When  third-order  con¬ 
stants  vanish  (i.e.,  a  St.  Venant-Kirchhoff  solid),  intrinsic 
stability  under  large  compression  tends  to  decrease  as 
Poisson’s  ratio  increases.  When  third-order  constants  are 
nonzero,  a  negative  (positive)  pressure  derivative  of  the 
shear  modulus  often  leads  to  unstable  modes  in  compres¬ 
sion  (tension).  When  pressure  derivatives  of  bulk  and 
shear  moduli  are  equal,  but  particular  third-order  elastic 
constants  differ,  internal  stability  is  the  same  for  spherical 
deformation  but  may  differ  for  uniaxial  tension/compres¬ 
sion.  For  simple  shear,  larger  magnitudes  of  third-order 
constants  promote  less  stable  behavior,  regardless  of  the 
sign  of  the  pressure  derivative  of  the  shear  modulus.  A 
compressible  neo-Flookean  elastic  model  demonstrates 
greater  intrinsic  stability  than  the  third-order  elastic  mod¬ 
el  when  Poisson’s  ratio  is  non-negative. 

Appendix  A.  Minimum  eigenvalue  ratios  in  the  limit  of 
incompressiblilty 

As  Poisson’s  ratio  v  — >  |,  an  isotropic  solid  approaches 
incompressibility.  Consider  a  Green  elastic  material  of  the 


type  described  in  Section  2.2,  with  vanishing  third-order 
elastic  constants.  At  a  fixed  strain  state,  define  ratios  of 
minimum  eigenvalues  as  Poisson’s  ratio  is  varied  by  the 
series 


^x|v=Q49  x 

AxUV  2 


Ax  |  y— 0499 
Ax  |  v-0.49 


Ax|v=0  5_iq-(*+1) 

Ax|v=o.5-Krk 

(A.1) 


Here,  X  can  refer  to  the  minimum  eigenvalue  of  stiffness 
matrix  A  or  B,  and  k  is  a  positive  integer.  Calculations  show 
that  r£  — ►  10  with  increasing  k  for  any  of  the  deformation 
paths  and  values  of  imposed  strain  e  considered  in 
Section  3,  and  for  both  X  =  A  and  X  =  B.  Representative 
results  are  listed  in  Table  3  for  spherical  and  uniaxial  com¬ 
pression  and  X  =  B,  demonstrating  how  rapidly  intrinsic 
stability  is  compromised  in  compression  as  the  material 
loses  compressiblity.  In  other  words,  since  <  0,  mini¬ 
mum  eigenvalues  become  strongly  negative  as  k  increases. 
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